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Abstract 

In the framework of relativistic SU(2)f baryon chiral perturbation theory we calcu- 
late the volume dependence of the nucleon mass up to and including 0(p 4 ). Since 
the parameters in the resulting finite size formulae are fixed from the pion mass de- 
pendence of the large volume nucleon masses and from phenomenology, we obtain 
a parameter-free prediction of the finite size effects. We present mass data from the 
recent Nf = 2 simulations of the UKQCD and QCDSF collaborations and compare 
these data as well as published mass values from the dynamical simulations of the 
CP-PACS and JLQCD collaborations with the theoretical expectations. Remark- 
able agreement between the lattice data and the predictions of chiral perturbation 
theory in a finite volume is found. 
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1 Introduction 



The computation of hadron masses is one of the basic goals of lattice QCD. 
However, as in the case of other observables, this computation suffers from a 
number of systematic uncertainties (not including the statistical errors): Lat- 
tice spacing and volume are finite, and the quark masses in the simulations are 
larger than in reality, in the extreme case of the sea quarks in the quenched 
approximation even infinite. Hence the Monte Carlo results must be supple- 
mented by several extrapolations: the continuum extrapolation a — ■> 0, the 
extrapolation to the thermodynamic limit L — > oo, and the chiral extrapola- 
tion sending the quark masses to their physical values. 

It is therefore an important issue to study hadron masses as functions of the 
quark masses or, in order to avoid problems connected with the evaluation of 
quark masses on the lattice, as functions of the pseudoscalar ("pion") mass. 
A suitable tool for investigating the quark-mass dependence of physical ob- 
servables like hadron masses is chiral perturbation theory (or, more generally, 
chiral effective field theory). Once one has convinced oneself that chiral per- 
turbation theory is applicable for the masses used in the simulations, one can 
extract low-energy constants of chiral effective field theory and extrapolate re- 
liably towards the physical masses. Indeed, as has recently been demonstrated, 
relativistic baryon chiral perturbation theory leads to a good chiral extrapo- 
lation function for the nucleon mass [1] (see also Ref. [2]) connecting available 
lattice results with the physical value. On the other hand, it is well known that 
finite size effects can be calculated from the same effective field theories that 
describe the quark-mass dependence [3,4]: As long as the volume is not too 
small, the finite size effects originate from pions which "propagate around the 
spatial box" . This approach leads to the so-called p expansion, valid for small 
pion masses m n = 0(p) and large spatial volumes L 3 with L~ x = 0(p) such 
that m n L = 0(p°). In this paper we work out the finite size corrections for 
the nucleon mass in this framework and present pion and nucleon masses ob- 
tained by the UKQCD and QCDSF collaborations in simulations with Nf = 2 
dynamical fermions. These and other recent nucleon mass data are compared 
with the formulae from chiral perturbation theory. Preliminary results of our 
investigation have been presented in Ref. [5]. For similar studies of the pion 
mass and pseudoscalar decay constants see Refs. [6,7]. Results from a different, 
but related, chiral analysis have been given in Ref. [8]. 



2 Finite size effects at 0(p 3 ) 

We follow Ref. [1] and employ relativistic SU(2)f baryon chiral perturbation 
theory as described in Ref. [9] . In this field theory the "elementary" degrees of 
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freedom are the pion and nucleoli fields. For the Lagrangian and further details 
see Refs. [10,1,9]. In particular, we shall use the so-called infrared regularisa- 
tion [9] , a variant of dimensional regularisation. The leading order contribution 
to the shift of the nucleon mass away from its value in the chiral limit comes 
from the piece in the 0{p 2 ) Lagrangian that breaks chiral symmetry explicitly. 
The next-to-leading order (NLO) contribution, i.e. the 0(p 3 ) contribution, is 
generated by the one-loop graph (a) of Fig. 1. 
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Fig. 1. One-loop graphs of NLO (a) and NNLO (b, c) contributing to the nucleon 
mass shift. The solid circle denotes a vertex from the leading order Lagrangian, the 
diamond a vertex from the 0(p 2 ) Lagrangian. 

The leading one-loop formula for the nucleon mass reads (in the infinite vol- 
ume) 



m N = m - Ac x m\ + 

Here and in the following, the constants g^, fn, ■ ■ ■ are to be taken in the chiral 
limit, mo denotes the nucleon mass in the chiral limit, and the quark mass 
has been replaced by the pion mass m^ using the Gell-Mann-Oakes-Renner 
relation. The pion decay constant f n is normalised such that its physical value 
is 92.4 MeV. The counterterm e^(A) is taken at the renormalisation scale A, 
which makes the result (1) scale independent. 

Expanding in powers of m^ (up to logarithms) we find 
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In Ref. [1] it was shown that this expansion is a good approximation of the 
full expression (1). Note, however, that all the terms oc m 3 , m*, ... in Eq. (2) 
are of the same chiral order. 
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At 0(p 3 ) the volume dependence of the nucleon mass comes from graph (a) 
in Fig. 1, where the pion couples to the nucleon through the pseudovector 
derivative coupling of the leading order Lagrangian [10]. The corresponding 
contribution to the nucleon mass reads in Minkowski space 



m a = -l- 



. 3g 2 A m ml 
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(3) 



Note that the rules of infrared regularisation have led us to take the integral 
over the Feynman parameter x from to oo rather than from to 1. After a 
Wick rotation this becomes in Euclidean notation 
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(5) 



In a finite spatial volume of linear size L (the time direction is left infinite) the 
integral over the spatial components of the loop momentum p is replaced by 
a sum over the discrete set of momenta allowed by the boundary conditions. 
Since we use periodic boundary conditions, the allowed momenta are of the 
form p = (2tt/L)£, where £ is a vector with integer components. Hence we get 
for the difference between the nucleon mass in a volume of size L 3 and the 
infinite volume nucleon mass at 0(p 3 ) 



m N (L) - m N (oo) = A a (L) = 



j—^ [p +p 4 + m x +m n (l 



x 



(6) 



According to Ref. [4] the difference between the sum and the integral is finite 
and given by 
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where the sum extends over all vectors ft with integer components excluding 
n = 0. Here rii can be interpreted as the number of times the pion crosses the 
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"boundary" of the lattice in the i direction. With the help of 

-oo / U \(l+o)/2 

j Q d\\ a e' XM2 - b / x = 2i^\ K 1+a (2VbAP) 
(see formula 3.471.9 in [11]) we obtain 
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where K v (x) is a modified Bessel function. For the numerical evaluation at 
smaller masses it is more advantageous to use the relation 
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For r = 2 we have from (9) 
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Using 

dx {x 2 + a 2 )- 1/4 K 1/2 {bVx^ + a 2 ) = J^-K (ab) (14) 
(see formula 6.596.3 in [11]) we obtain finally 

m N (L) - m N (oo) = A (L) 
Sg^moml 



16tT 2 / 2 



^°°cb X^'i^o (L|n|^/mgx 2 + m2(l -a;)) . (15) 



This is the complete 0(p 3 ) result for the volume dependence of the nucleon 
mass in relativistic baryon chiral perturbation theory. Note that the finite 
volume has not introduced any new parameter. 



3 Finite size effects at 0(p 4 ) 



At 0(p 4 ), the nucleon self-energy receives additional contributions from graphs 
(b) and (c) in Fig. 1. Up to higher-order corrections, the contribution to the 
nucleon mass from graph (c) is cancelled by the contribution arising from the 
insertion of = m — 4cim 2 in the 0(p 3 ) piece corresponding to graph 
(a). Alternatively, one could replace the mass m in the free Lagrangian by 
m — 4cim 2 and omit graph (c) [9]. Thus only graph (b) is relevant for us. The 
resulting formula for the nucleon mass in the infinite volume at 0(p 4 ) involves 
two new coupling constants c 2 and c 3 and reads [1,9] 
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The contribution to the nucleon mass from graph (b) in Fig. 1 is given by 
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in Minkowski space. After a Wick rotation this becomes in Euclidean notation 
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Hence we get an additional contribution to the difference between the nucleon 
mass in a volume of size L 3 and the infinite volume nucleon mass: 
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Now we can use Eq. (9) to rewrite this result as 
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and by means of Eq. 6.596.3 in [11] we get the final formula 
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Hence we have 



m N (L) - m N (oo) = A a (L) + A b (L) + 0(p 5 ) , (24) 

where A (L) is given in Eq. (15). Eq. (24) represents the complete 0(p 4 ) result 
for the volume dependence of the nucleon mass in relativistic baryon chiral per- 
turbation theory and will serve as the basis for our numerical studies. Again, 
the finite volume did not lead to the appearance of additional parameters. 

The leading contribution in Eq. (24) comes from the terms with \n\ = 1 
corresponding to pions which travel around the spatial box exactly once. In the 
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Appendix we bring this contribution into the "dispersive" form introduced by 
Liischer [12], which relates the finite volume effects to the pion-proton forward 
elastic scattering amplitude. We find 
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The coefficients 2ci — c 3 and c 2 are related to coefficients in the so-called 
subthreshold expansion [13]. Unless m^L is rather large, the subleading terms 
with \n\ > 1 are not negligible. 



4 Monte Carlo data 

We want to compare our formulae to Monte Carlo data for nucleon masses. 
As the finite size effects rely essentially on the existence of a "pion cloud" 
and hence on the presence of sea quarks, such a comparison makes sense only 
for masses extracted from dynamical simulations. Large scale simulations with 
two flavours of dynamical quarks are beginning to deliver results, but the quark 
masses are still rather large, and the continuum limit is not easy to perform 
reliably, because only a rather small range in the lattice spacing a is covered. 
Therefore it is not straightforward to set the scale for the simulation results. A 
popular method is to use the force parameter r [14] for that purpose, because 
it is expected to depend only weakly on the quark masses. In the following we 
shall adopt this procedure using r = 0.5 fm. However, it must be noted that 
recently some doubts have been raised on the reliability of this method [15], 
and it would certainly be advantageous to avoid the scale problem completely 
by considering dimensionless ratios, e.g., ratios of masses with the pion decay 
constant and performing the whole analysis for these quantities. Another point 
of concern is the question of whether the quark masses in the simulations are 
small enough to justify the application of chiral perturbation theory. 

Over the last years, the UKQCD and QCDSF collaborations have generated 
gauge field configurations with two flavours of dynamical quarks for a variety 
of quark masses, lattice spacings and volumes. These simulations use the stan- 
dard Wilson plaquette action for the gauge fields and the non-perturbatively 
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0(a) improved clover action for the fermions. The algorithm employed is the 
Hybrid Monte Carlo algorithm, which recently could be sped up consider- 
ably [16]. Details of the extraction of r /a are given in [17]. 

The computation of the nucleon masses on the configurations generated by 
UKQCD is also described in Ref. [17]. On the configurations generated by 
QCDSF a somewhat different procedure has been used. While standard in- 
terpolating fields have been employed in both cases, QCDSF has applied Ja- 
cobi smearing instead of fuzzing when computing the quark propagators from 
point sources. The nucleon masses have been obtained by fitting the source and 
sink smeared correlation function to exp [— m^t] + A^* exp [— m^* (T — t)] , 
where N* is the negative parity partner of the nucleon. The x 2 was calculated 
from the diagonal part of the covariance matrix only. But it has been checked 
that using the full covariance matrix gives consistent results. The fit range 
was fixed by looking for a region where the results for am N and am N * are 
independent of the fit range. 

Table 1 

Simulation parameters and results from the UKQCD and QCDSF collaborations. 





Coll. 







^sea 


volume 


r /a 




arriN 


1 


QCDSF 


5.20 


0. 


1342 


16 3 


X 


32 


4.077(31) 


0.5841(11) 


1.1071(40) 


2 


UKQCD 


5.20 


0. 


1350 


16 3 


X 


32 


4.754(40) 


0.405(5) 


0.883(12) 


3 


UKQCD 


5.20 


0. 


1355 


16 3 


X 


32 


5.041(40) 


0.294(4) 


0.766(11) 


4 


UKQCD 


5.20 


0. 


13565 


16 3 


X 


32 


5.246(51) 


0.2470(40) 


0.676(24) 


5 


UKQCD 


5.20 


0. 


1358 


16 3 


X 


32 


5.320(50) 


0.2080(70) 


0.636(33) 


6 


QCDSF 


5.25 


0. 


1346 


16 3 


X 


32 


4.737(21) 


0.4925(16) 


0.9455(57) 


7 


UKQCD 


5.25 


0. 


1352 


16 3 


X 


32 


5.138(45) 


0.3842(16) 


0.8032(69) 


8 


QCDSF 


5.25 


0. 


13575 


24 3 


X 


48 


5.430(60) 


0.2599(15) 


0.6160(58) 


9 


UKQCD 


5.26 


0. 


1345 


16 3 


X 


32 


4.708(52) 


0.509(2) 


1.011(10) 


10 


UKQCD 


5.29 


0. 


1340 


16 3 


X 


32 


4.813(45) 


0.577(2) 


1.086(9) 


11 


QCDSF 


5.29 


0. 


1350 


16 3 


X 


32 


5.227(37) 


0.4208(8) 


0.8344(34) 


12 


QCDSF 


5.29 


0. 


1355 


12 3 


X 


32 


5.756(33) 


0.3637(48) 


0.864(12) 


13 


QCDSF 


5.29 


0. 


1355 


16 3 


X 


32 


5.560(30) 


0.3334(15) 


0.7188(76) 


14 


QCDSF 


5.29 


0. 


1355 


24 3 


X 


48 


5.566(20) 


0.3265(6) 


0.6857(35) 



In addition we use data available in the literature [18,19]. The CP-PACS col- 
laboration [18] works with a renormalisation-group improved gauge action 
and a mean field improved clover quark action, while the JLQCD collabora- 
tion [19] employs the same actions as UKQCD and QCDSF. However, JLQCD 
uses a slightly different value for the improvement parameter cgw- In Table 1 
we present our results including those that have already been published in 
Ref. [17] (points 2, 3, 9, 10). Relevant data obtained by the CP-PACS and 
JLQCD collaborations are collected in Table 2. Some important quantities 
converted to physical units are given in Table 3. 
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Table 2 

Simulation parameters and results from the CP-PACS [18] and JLQCD [19] collab- 
orations. 





Coll. 





^sca 


volume 


Tc\l a 






15 


CP-PACS 


1.95 


o. 


1410 


16 3 


x 


32 


3.014(33) 


0.42700(98) 


1.0532(51) 


16 


CP-PACS 


1.95 


n 


1 400 


16 3 


x 


32 


2.821(29) 


0.59580(69) 


1.2679(39) 


17 


CP-PACS 


1.95 


n 


1 390 


16 3 


x 


32 


2 651(42) 

___J • V ' ~.S ____. \ A ___J 1 


72857(68) 


1 4559(38) 


18 


CP-PACS 


1.95 





1375 


16 3 


x 


32 


2.497(54) 


0.89400(52) 


1.7035(34) 


19 


CP-PACS 


2.10 





1382 


24 3 


x 


48 


4.485(12) 


0.29459(85) 


0.7204(42) 


20 


CP-PACS 


2.10 


0. 


1374 


24 3 


x 


48 


4.236(14) 


0.42401(46) 


0.8955(35) 


21 


CP-PACS 


2.10 


o. 


4367 


24 3 


x 


48 


4.072(15) 


0.51671(67) 


1.0226(32) 


22 


CP-PACS 


2.10 


0. 


1357 


24 3 


x 


48 


3.843(16) 


0.63010(61) 


1.1855(26) 


23 


CP-PACS 


2.20 


0. 


1368 


24 3 


x 


48 


5.410(21) 


0.2785(22) 


0.6314(55) 


24 


CP-PACS 


2.20 


0. 


1363 


24 3 


x 


48 


5.237(22) 


0.3554(10) 


0.7349(42) 


25 


CP-PACS 


2.20 


0. 


1358 


24 3 


X 


48 


5.073 19 


0.4190(13) 


0.8252(47) 


26 


CP-PACS 


2.20 


0. 


1351 


24 3 


X 


48 


4.913(21) 


0.49996(83) 


0.9330(76) 


27 


JLQCD 


5.20 


0. 


1340 


12 3 


X 


48 


3.826(50) 


0.619(10) 


1.153(12) 


28 


JLQCD 


5.20 


0. 


1343 


12 3 


X 


48 


4.031(84) 


0.5474(51) 


1.094(23) 


29 


JLQCD 


5.20 


0. 


1346 


12 3 


X 


48 


4.200(56) 


0.5011(70) 


1.006(24) 


30 


JLQCD 


5.20 


0. 


1350 


12 3 


X 


48 


4.481(67) 


0.4239(71) 


0.915(34) 


31 


JLQCD 


5.20 


0. 


1355 


12 3 


X 


48 


5.06(12) 


0.328(14) 


0.820(36) 


32 


JLQCD 


5.20 


0. 


1340 


16 3 


X 


48 


3.880(36) 


0.6200(21) 


1.166(20) 


33 


JLQCD 


5.20 


0. 


1343 


16 3 


X 


48 


4.098(45) 


0.5528(40) 


1.063(13) 


34 


JLQCD 


5.20 


0. 


1346 


16 3 


X 


48 


4.287(58) 


0.4939(20) 


0.991(14) 


35 


JLQCD 


5.20 


0. 


1350 


16 3 


X 


48 


4.621(42) 


0.4003(33) 


0.828(12) 


36 


JLQCD 


5.20 


0. 


1355 


16 3 


X 


48 


5.059(71) 


0.2806(64) 


0.707(29) 


37 


JLQCD 


5.20 


0. 


1340 


20 3 


X 


48 


3.946(30) 


0.61630(55) 


1.1566(26) 


38 


JLQCD 


5.20 


0. 


1343 


20 3 


X 


48 


4.143(29) 


0.55270(62) 


1.0626(23) 


39 


JLQCD 


5.20 


0. 


1346 


20 3 


X 


48 


4.336(50) 


0.49020(71) 


0.9644(31) 


40 


JLQCD 


5.20 


0. 


1350 


20 3 


X 


48 


4.635(53) 


0.40037(55) 


0.8252(26) 


41 


JLQCD 


5.20 


0. 


1355 


20 3 


X 


48 


5.092(83) 


0.27133(72) 


0.6468(36) 
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Table 3 

Simulation data converted to physical units using vq = 0.5fm. 





a [fm] 


L [fm] 




m n [GeV] 


m N [GeV] 


1 


0.12 


1.96 


9.3 


0.9398(74) 


1.781(15) 


2 


0.11 


1.68 


6.5 


0.760(11) 


1.657(26) 


3 


0.10 


1.59 


4.7 


0.5849(92) 


1.524(25) 


4 


0.10 


1.52 


4.0 


0.5114(97) 


1.400(52) 


5 


0.09 


1.50 


3.3 


0.437(15) 


1.335(70) 


6 


0.11 


1.69 


7.9 


0.9207(51) 


1.768(13) 


7 


0.10 


1.56 


6.1 


0.7791(76) 


1.629(20) 


8 


0.09 


2.21 


6.2 


0.5570(69) 


1.320(19) 


9 


0.11 


1.70 


8.1 


0.946(11) 


1.878(28) 


10 


0.10 


1.66 


9.2 


1.096(11) 


2.063(26) 


11 


0.10 


1.53 


6.7 


0.8681(64) 


1.721(14) 


12 


0.09 


1.04 


4.4 


0.826(12) 


1.963(29) 


13 


0.09 


1.44 


5.3 


0.7316(51) 


1.577(19) 


14 


0.09 


2.16 


7.8 


0.7172(29) 


1.5062(94) 


15 


0.17 


2.65 


6.8 


0.5079(57) 


1.253(15) 


16 


0.18 


2.84 


9.5 


0.6633(69) 


1.412(15) 


17 


0.19 


3.02 


11.7 


0.762(12) 


1.523(24) 


18 


0.20 


3.20 


14.3 


0.881(19) 


1.679(36) 


19 


0.11 


2.68 


7.1 


0.5214(21) 


1.2751(82) 


20 


0.12 


2.83 


10.2 


0.7088(25) 


1.4971(77) 


21 


0.12 


2.95 


12.4 


0.8304(32) 


1.6434(79) 


22 


0.13 


3.12 


15.1 


0.9556(41) 


1.7980(85) 


23 


0.09 


2.22 


6.7 


0.5946(52) 


1.348(13) 


24 


0.10 


2.29 


8.5 


0.7345(37) 


1.519(11) 


25 


0.10 


2.37 


10.1 


0.8389(41) 


1.652(11) 


26 


0.10 


2.44 


12.0 


0.9694(45) 


1.809(17) 


27 


0.13 


1.57 


7.4 


0.935(19) 


1.741(29) 


28 


0.12 


1.49 


6.6 


0.871(20) 


1.740(52) 


29 


0.12 


1.43 


6.0 


0.831(16) 


1.667(46) 


30 


0.11 


1.34 


5.1 


0.750(17) 


1.618(65) 


31 


0.10 


1.19 


3.9 


0.655(32) 


1.637(82) 


32 


0.13 


2.06 


9.9 


0.9494(94) 


1.785(35) 


33 


0.12 


1.95 


8.8 


0.894(12) 


1.719(28) 


34 


0.12 


1.87 


7.9 


0.836(12) 


1.677(33) 


35 


0.11 


1.73 


6.4 


0.7300(90) 


1.510(26) 


36 


0.10 


1.58 


4.5 


0.560(15) 


1.412(61) 


37 


0.13 


2.53 


12.3 


0.9598(74) 


1.801(14) 


38 


0.12 


2.41 


11.1 


0.9037(64) 


1.737(13) 


39 


0.12 


2.31 


9.8 


0.8388(98) 


1.650(20) 


40 


0.11 


2.16 


8.0 


0.7324(84) 


1.509(18) 


41 


0.10 


1.96 


5.4 


0.5453(90) 


1.300(22) 
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5 Comparison with chiral perturbation theory 

In this section we confront the nucleon masses presented in the preceding 
section with chiral perturbation theory. We use data obtained on (relatively) 
large lattices to find appropriate values for the parameters of the chiral expan- 
sion. After fixing these parameters, the finite size formulae predict the volume 
dependence of the nucleon masses without any free parameter. 

In Ref. [1] an analysis of data points from Tables 1 and 2 obtained on (rela- 
tively) large and fine lattices has been presented (see also Ref. [2]). More pre- 
cisely, masses from simulations with a < 0.15 fm, m^L > 5 and m n < 800 MeV 
have been considered. Here we shall impose the same conditions and select the 
data points 2, 7, 8, 14, 19, 20, 23, 24, 40, 41 for the further analysis. (Point 13, 
which also fulfills the above criteria, is discarded, because the corresponding 
simulation parameters agree with those of point 14, except that the lattice is 
smaller. Similarly, simulations 30 and 35 just repeat simulation 40 on smaller 
lattices.) 

We cannot determine all parameters in the chiral perturbation theory formu- 
lae from fits to the lattice data. Therefore we decided to fix g^, fn, C2, C3, and 
to fit the remaining parameters uiq, ci, e{(\) choosing A = 1 GeV. For qa and 
f n we take the physical values qa = 1.267, /„- = 92.4 MeV. Choosing values 
for C2 and C3 is more difficult, although there are quite a few phenomeno- 
logical determinations in the literature (see, e.g., [20,21,22,23,24,25,26,27,28]). 
In particular, C3 seems to be subject to a considerable uncertainty. We set 
C2 = 3.2 GeV -1 , which is compatible with Refs. [20,21,22,23,27], and consider 
two possibilities for C3: C3 = —3.4 GeV" 1 and C3 = — 4.7 GeV" 1 . According 
to Refs. [27,28], the value c 3 = —3.4 GeV -1 is consistent with the empirical 
nucleon-nucleon phase shifts and with the value extracted in [24] from pion- 
nucleon scattering. On the other hand, C3 = —4.7 GeV" 1 is the central value 
obtained in the pion- nucleon scattering analysis of Ref. [24] , albeit with large 
error bars, and is consistent with the results of Ref. [26] . 

From a fit of Eq. (16) to the ten points selected above with c 3 = —3.4 GeV -1 
we obtain the results labelled as "Fit 1" in Table 4. Data points and fit are 
shown in Fig. 2. (Using qa = 1.2, f n = 88 MeV as better approximations 
to the values in the chiral limit does not lead to significant differences.) It 
is first of all remarkable that the Monte Carlo data obtained by different 
collaborations with different actions all lie close to a single curve. Thus there 
seems to be no sign of appreciable lattice artefacts. Secondly, the fit describes 
the data very well and is at the same time compatible with the physical mass 
values. Thirdly, the fit parameters take values which are consistent with known 
phenomenology. In particular, c\ compares favourably with phenomenological 
determinations [20,23,24,25]. Repeating the fit with c 3 = —4.7 GeV" 1 leads 
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to the results labelled as "Fit 2" in Table 4, which are farther away from the 
phenomenological numbers. 

Table 4 

Results from the four fits of nucleoli mass data as described in the text. The pa- 
rameter C3 has been kept fixed. 





c 3 [GeV- 1 ] 


m [GeV] 


ci [GeV- 1 ] 


el (A = lGeV) [GeV" 3 ] 


x 2 


Fit 1 


-3.4 


0.89(6) 


-0.93(5) 


2.8(4) 


12.18 


Fit 2 


-4.7 


0.76(6) 


-1.25(5) 


1.7(5) 


11.85 


Fit 3 


-3.4 


0.88(-) 


-0.93(4) 


3.0(6) 


0.29 


Fit 4 


-4.7 


0.87(-) 


-1.11(4) 


3.2(6) 


0.39 




Fig. 2. Nucleon mass data on (relatively) large and fine lattices. The star indicates 
the physical point. The curve corresponds to Fit 1 in Table 4. 

Of course, one may doubt the applicability of chiral perturbation theory at the 
larger pion masses which enter our analysis. On the other hand, we do not have 
enough results at smaller masses to restrict ourselves to a safer mass range 
performing the same kind of fits as above. We can, however, constrain the fit 
by requiring the fit curve to pass through the physical point. We implement 
this constraint by choosing m such that the condition is satisfied. In this way 
we eliminate one fit parameter leaving only c\ and e{(\ = 1 GeV) to be fitted. 
The reduced number of parameters then allows us to restrict the fit to a smaller 
number of masses. From such constrained fits to the four data points at the 
smallest masses (points 8, 19, 23, 41) we get the results labelled as "Fit 3" 
(with c 3 = -3.4 GeV -1 ) and "Fit 4" (with c 3 = -4.7 GeV" 1 ) in Table 4. While 



13 



the results of Fit 1 and Fit 3 are well compatible with each other (and with 
phenomenology), the constraint has a stronger effect for c 3 = — 4.7GeV~ 1 . 

Using the parameters obtained from Fit 1, which are consistent with phe- 
nomenology and the large volume Monte Carlo data, we can now evaluate our 
finite size corrections. Note that no further free parameters are involved. For 
the comparison with the Monte Carlo results we need data from simulations 
which differ only in the lattice size while all other parameters are kept fixed. 
There are six sets of three simulations each among the results in Tables 1, 2. 
We choose the three sets with the smallest pion masses and plot the nucleon 
masses versus the lattice size in Figs. 3, 4, 5. The curves are computed from 
our finite size formulae. The curve labelled p 4 corresponds to 

m N (L)=A a (L) + A b (L)+m N (oo), (26) 

where m^oo) is determined such that 7Bat(L) on the largest lattice agrees 
with the Monte Carlo value. For the pion mass we use the value from the 
largest lattice. In the curve labelled p 3 , the p 4 contribution A b (L) has been 
left out, but mjv(oo) has not been changed. 



2.2 




3 

P 4 


2.0 




P 


1.8 






£ 1-6 
O 






|1.4 






1.2 






1.0 






0.8 


1,1 





3 I I I I I I I I I I I 

0.0 0.5 1.0 1.5 2.0 2.5 



L [fm] 

Fig. 3. Volume dependence of the nucleon mass for m n = 545 MeV (data points 31, 
36, 41). The dotted curve shows the contribution of the p 3 term, while the solid 
curve includes also the p 4 correction, with the parameters taken from Fit 1. 

Our formula describes the volume dependence of the nucleon mass remarkably 
well. However, the agreement would deteriorate had we taken into account only 
pions which travel around the lattice exactly once, by restricting the sums over 
n in Eqs. (15), (23) to n with \n\ = 1. 
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2.2 
2.0 
1.8 



% 1-6 
O 

£ 1.4 



1.2 
1.0 



0.8 




0.0 0.5 1.0 1.5 2.0 2.5 

L [fm] 



Fig. 4. Volume dependence of the nucleon mass for = 717 MeV (data points 12, 
13, 14). The dotted curve shows the contribution of the p 3 term, while the solid 
curve includes also the p 4 correction, with the parameters taken from Fit 1. 



2.2 
2.0 
1.8 
% 1-6 

o 

I 1.4 
1.2 
1.0 



0.8 




0.0 0.5 1.0 1.5 2.0 2.5 

L [fm] 



Fig. 5. Volume dependence of the nucleon mass for = 732 MeV (data points 30, 
35, 40). The dotted curve shows the contribution of the p 3 term, while the solid 
curve includes also the p 4 correction, with the parameters taken from Fit 1. 
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6 Conclusions 



Chiral effective field theories evaluated in a finite volume yield expressions for 
physical quantities, e.g. masses, in the form of an expansion with an expansion 
parameter p, where the pion mass m n and the inverse box size XjL both count 
as quantities of order p. We have applied this formalism to the nucleon mass 
using relativistic SU(2)f baryon chiral perturbation theory up to and including 
0(p 4 ) terms. 

Setting the scale with r = 0.5 fm, we can convert Monte Carlo results obtained 
in simulations with dynamical quarks to physical units and compare them 
with our formulae. The infinite-volume formulae lead to a surprisingly good 
description of the data coming from (relatively) large lattices, with values for 
the coupling constants of the effective theory which are completely consistent 
with known phenomenology [1,2]. Given these coupling constants, there are no 
free parameters left in the finite size formulae. But the p A formulae reproduce 
the finite size effects observed in the Monte Carlo data surprisingly well. For 
this agreement it is essential that our formalism incorporates the effects of 
pions travelling around the lattice arbitrarily many times. Thus we arrive at 
a more optimistic opinion on the applicability of chiral perturbation theory 
than previous investigations [29]. 

We find a remarkably consistent picture of the volume and pion mass depen- 
dence of the nucleon mass, in spite of the fact that the "convergence" of the 
chiral expansion at the rather large pion masses used in the available sim- 
ulations is far from obvious. Lattice data for smaller masses will therefore 
be of great importance to corroborate this conclusion. Another point which 
deserves further investigation is the cutoff and mass dependence of r [15]. 
Alternatively one could try to avoid the use of r by basing the analysis on 
ratios of masses with a quantity which is more easily accessible to chiral per- 
turbation theory than r , e.g. the pion decay constant. One could then work 
out formulae directly for these ratios. Finally it would be interesting to see if 
the pion mass can be described along similar lines. 



Appendix 

In this Appendix we compare our description of finite size effects with a finite 
size formula derived by Liischer [12]. He expresses the leading contributions 
to the relative finite size effect 




(27) 
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in terms of the pion-nucleon coupling constant g^N with the physical value 

^ « 14.3 (28) 

and the pion-proton forward elastic scattering amplitude F np (v). Here v is the 
"crossing variable" v = (s — u)/(Am N ) with m N = m N (oo). Liischer finds 



- sdb ©71^ exp (--^VW) ^(i^P) 

+ O (V am - L ) , (29) 



where a > y3/2. The amplitude F np {v) is related to the amplitude C + (z/) = 
of Ref. [13] through 

= 6m w C + (z/). (30) 

It can conveniently be decomposed into the pseudovector Born term and a 
remainder R(v): 



6g 



F *pM = i r4^T7 4 + Riy) ■ ( 31 ) 

1 — Am N v 2 1 m\ 

Using the so-called subthreshold expansion [13] for vanishing momentum trans- 
fer 

oo 

R(u) = 6m N J2dt u 2k (32) 

k=0 

Liischer's formula becomes (ignoring questions of convergence) 



rLuscher 9 { m w \ 2 2 \ 71 ( ml 

$N = 7TT &n{ o 7 eX P - m nL\ l ~ 



Q7r z Vmjv/ Im^L \ V 4m^ 



1 /-oo exp ( — m^Ly/l + p 2 
dp 



m n L J-oo r l + Am 2 N p 2 /m 2 
9 rn„ ^__ 2fcj+ , lNfe /-dp 



+ O (V am - L ) . (33) 

As it is written, Liischer's formula refers to a fixed value of m n . However, as 
far as we have information on the quark- mass dependence of mN, g n N and the 
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scattering amplitude - and chiral perturbation theory provides such informa- 
tion - we can make use of this knowledge to describe the m n dependence of 
<5jv as well. But Liischer's formula takes into account only pions which travel 
around the lattice exactly once. For the comparison with our formulae we 
therefore have to restrict the sums over n in Eqs. (15), (23) to n with \n\ = 1. 
In order to get the relative finite size effect we must furthermore divide by 
Tfijy . Because our finite size formulae correspond to 0(p 3 ) (in the case of A ) 
and 0(p 4 ) (if also A& is taken into account) in the chiral expansion, we identify 
here and in the following m N with m - the higher terms in mjv would give a 
contribution only at 0(p 5 ). Thus we get from graph (a) 



= 3g|m| 
° 167T 2 /. 2 

while graph (b) yields 



J dx ■ 6 • Kq ^L\JmQX 2 + rn%( 



1 — x~ 



(34) 



S h = 



ml 



2Tr 2 f 2 m 



Ki(m n L) K 2 (m n L) 
(2ci - c 3 ) h c 2 



(35) 



First we consider 5 a and show that the integral over x in (34) can be rewritten 

as 



J dx Kq ^L^jmQX 2 + m 2 (l — x) 



7T 



m L 



exp —m-xLi 1 — 



dp 



exp (—m n Ly/l +p 2 ) 



4ml J m n L J ~oo 1 + 4m\ p 2 / 'm\ 



(36) 



for m\ < 4m 2 ,. We begin by writing 



J dx Kq ^L\Jmlx 2 + m 2 (l — x) 
= j dx Kq (^L\JmQX 2 + m 2 (l — x) 

— J dx Kq [L\jm\x 2 + m 2 (l — ~x)j . 



(37) 



The first integral on the right-hand side can be evaluated with the help of 
formula 6.596.3 in [11]: 
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oo 

dxK I L A 

OO 



* ( t L m l 

exp —m n L\l — - — k 



m§L \ "V Ami J 
In the second integral we use the integral representation 

/oo , . — 1/2 

ds (s 2 - y 2 ) e~ xs 

(valid for x, y > 0) and interchange the order of the integrations: 



( dxK (L^mlx 2 + ml(l — x) 

J —oo \ 



Jl Jo \ rn2 



mi 



with the abbreviation 



V[S) moV + 4m 2 2mg - 

The integration over x can now be performed leading to 

'2m n 



^ fdse-— arctanf^v^T) . 
m Ji V to^ / 



Partial integration yields 



2 /■oo p— ra^Ls 

ds 



and the change of variables s = \/l + p 2 gives the final expression 

rO / , \ 1 roo exp (—m w Ly/l + p 2 J 

/ dxiroUV^o^ + ^l 1 -^) = r/ d P — ; | , 2 2i 2 

■/-oo V v / m n L J -oo l + 4moP J /m; 

completing the proof of Eq. (36). 

With the help of the Goldberger-Treiman relation 

Pa — 9ttN — , 
m 

which is exact in the chiral limit, we now obtain 
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1 /-oo exp (—m n Ly/l + p 2 ) ) 

r / dp } - 2 — - 7 I. (46) 

m w L J -oo l + 4m 2 ,p 2 /m 2 J 

The chiral corrections to (45) are of order p 2 and hence would contribute only 
at 0(p 5 ), beyond the accuracy of our calculation. Thus <5 a reproduces the first 
two terms in (33) up to a factor of two in the first term, which in Liischer's 
treatment originates from nucleon exchange diagrams. 

Do the following terms of Eq. (33) contain our 5b, Eq. (35)? With the help of 
and the tree-level relations (see, e.g., Ref. [9]) 



2m 2 

4, = --7r(2ci-c 3 ) + c>("4), 

J 7T 

^0 = ^02 + 0(77^) (48) 



we obtain 



9 m„- 



9 mi 



27r 2 / 2 m 
in complete agreement with (35). 



K^m^L) K 2 (m 7r L) 
(2ci - c 3 ) h c 2 — rr^- + 



(49) 



In short, when we write the \n\ = 1 contributions of our finite size formula 
(24) in the form that Liischer uses [12], we find, within the accuracy of our 
chiral expansion, 



- sdb OVIS exp (-^vW) f ^p) > ( 5 °) 

which looks like Liischer's formula, except that the leading term has a coeffi- 
cient twice as large. Numerically this makes a significant difference (more than 
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a factor 2) because in Liischer's case the exponential and the integral almost 
cancel each other. 
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